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LONG  RANGE  TRAJECTORY  PREDICTION  ERRORS 
FOR  LEAST  SQUARES  SMOOTHING 


ABSTRACT 


This  paper  presents  exact  formulas  and  graphs  to 
determine  the  root  mean-square  error  of  prediction 
for  position  and  velocity  estimates  of  a  target, 
based  on  least-squares  smoothing.  The  prediction 
interval  is  taken  from  an  end-point  estimate  up  to 
ten  times  the  smoothing  interval.  The  r.m.s.  values 
may  be  used  to  determine  how  large  a  search  volume 
must  be  to  accommodate  a  substantial  data  dropout 
interval . 


29  September  1969 


3 


SP-3424 


1.0  INTRODUCTION 

The  application  considered  is  as  follows:  Suppose  one  observes 
tracking  data  over  M  sample  points,  spaced  T  seconds  apart,  and  smooths  the 
data  by  curve  fitting  a  least-squares  polynomial  of  degree  n  <  M  to  the  data. 

One  wishes  a  "best"  estimate  of  the  position  and  velocity  of  the  target  at  a 
time  TM(1  +  p)  as  measured  from  the  beginning  of  the  data  set.  Values  of  p 
considered  in  this  paper  are  shown  in  Figures  1  through  5  (See  Appendix)  for 
p  =  0,  .5  and  1  to  10.  Thus  if  tracking  data  is  obtained  over  say  M  =  10  samples 
for  T  =  1.  second  sampling  interval,  and  a  data  dropout  occurs  for  ten  seconds, 
one  could  curve  fit  a  least-squares  quadratic  or  cubic  to  the  ten  data  points, 
and  estimate  the  target's  position  and  velocity  ten  seconds  later.  The  r.m.s. 
of  the  estimate  can  be  used  to  determine  the  size  of  the  search  volume  to 
find  the  target. 

The  r.m.s.  values  must  be  considered  as  lower  bounds  to  determine 
search  volumes  since  bias  errors,  due  to  departure  of  the  actual  trajectory 
from  the  low-degree  polynomical,  will  increase  the  prediction  error. 

The  curves  presented  in  this  paper  for  the  r.m.s.  values  are  deter¬ 
mined  from  exact  equations.  A  set  of  asymptotic  values  of  r.m.s.  for 
M  >  >  1  and  p  >  >  1  are  also  derived, 

1.1  SOLUTION  FOR  MEAN  SQUARE  ERROR  (.5 2) 

Formulas  for  computing  the  r.m.s.  for  the  output  of  a  finite  memory 
optimal-digital  filter,  when  the  input  signal  consists  of  a  nonrandom  poly¬ 
nomial  of  degree  n  plus  uncorrelated  stationary  noise,  have  been  published  [1], 
The  optimal  filter  is  equivalent  to  least-squares  curve  fitting  of  an  n^ 
degree  polynomial. 
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It  is  shown  that: 


n  5 (K)  (M  +  a)  2 

Y  _ =}  5  6  (M,  a,  K,  n) 

L**K  S  (L,M) 


S  (L,M) 


AT  2  2 

(L!)  M  7T  (M  -  j  ) 


(2L) !  (2L  +  1)! 


(1) 


(2) 


where 

interval  between  samples, 
mean  square  of  input  noise  N(u). 

number  of  samples  in  digital  filter  of  fixed  memory  length 
(M  -  1)  T. 

mean  square  of  output  of  filter. 

derivative  of  desired  output, 
time  at  which  desired  output  is  evaluated, 
th 

K  derivative  of  the  orthogonal  polynomial  of  degree  L 
evaluated  at  t  =  (M  +  a)*T 

In  the  referenced  paper,  extensive  curves  6  are  presented  for  a  =  0, 

or  end-point  prediction  a  =  1,  or  next-point  prediction  and  a  =  -[M  -  l]/2, 

th  st  r»d 

or  smoothing  of  the  mid-point  of  the  data  interval  for  o  ,  1  ,  and  2 

derivative  estimates.  The  computations  for  <5  are  extended  for  long-range 
prediction  errors  in  the  following. 


T  = 

2 

°N  = 
M  = 

2 

°A  = 
K  = 

T*  (M  +  a)  = 
^  (M  +  a)  = 
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Let  the  estimate  be  desired  at  the  normalized  (T=.l)  time 
U  =  (1  +  p)  •  M 


corresponding  to  a  =  Mp.  Curves  of  6  vs  M  are  presp^t^d  for  values  of 
0  <_  p  _<  10,  n  =  1  to  3  for  zero  derivative  estimates  and,  n  =  2,  3,  for  first 
derivative  estimates.  It  is  easily  shown  that 


62(M,  Mp,  0,  0) 


1  =  1 
M  S (0 ,  M) 


(3) 


independent  of  p,  and 


62(M,  Mp,  1,  0)  =  0 


6  2  (M , 


Mo,  1,  1)  = 


12 


M(M 


1) 


S (L,  M) 


(4) 


independent  of  ,  . 

Asymptotic  formulas  are  presented  for  S  as  M  -*•  00  and  an  approximation  to 
6  are  given  for  moderate  values  of  M. 

1.2  APPROXIMATE  AND  ASYMPTOTIC  FORMULAS  FOR  5 2 

2 

In  this  section,  simplified  equations  for  5  will  be  derived  for 
K  =  0,  and  asymptotic  formulas  for  K  =  0  and  K  =  1. 


The  term  S(L,  M)  can  be  written 
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for  M  >>  1  n  >  1 


62  (M,  Mp,  0,  0)  -  1/M 


6 


2 


(M,  Mp,  0,  1)  S  1/M 


r 2  n 

h  +  l 


2 

6  (M,  Mp  ,  0,  2)  =1/M 


2 

6 


(M,  Mp,  0,  3)  S  1/M 


"  6 

n_ 

c3 


660n| 


(11) 


As  n,  M  -*•  <*  one  has  the  asymptotic  equations 


6  (M,  Mp,  0,  L)~  n  L  -  0,  1,  2,  ... 


(12) 


Similarly  for  first  derivative  estimation,  the  orthogonal  polynomials 
evaluated  at  U  are  given  by 


'1=  1 


f,2  =  2 [U  -  —-1] 


S  =  3[U-  2 


M  +  1,2  r3M  -  7 


--]  ~  [ 


20 


]. 


(13) 


For  M  >>  1,  n  i  1,  an  approximate  formula  for  <S  is  given  by 
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6  (M,  Mp,  1,  0)  =  0 


2  '12 
6  (M.  Mp,  1,  1)  =  — 

M 


S2(M,  Mp,  1,  2)  =  [1  +  60n2] 

MJ 

62(M,  Mp,  1,  3)  =  1-8-°~-  [I4r|2  -  1]. 

M 


For  M  >>  1,  n  >>  1  an  asymptotic  formula  for  6  is  given  by 


6 (M,  Mp,  1,  L) 
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